Cooper-Pair Injection into Quantum Spin Hall Insulators 
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We theoretically study tunneling of Cooper pairs from a superconductor spanning a two- 
dimensional topological insulator strip into its helical edge states. The coherent low-energy electron- 
pair tunneling sets off positive current cross correlations along the edges, which reflect an interplay 
of two quantum-entanglement processes. Most importantly, superconducting spin pairing dictates 
a Cooper pair partitioning into the helical edge liquids, which transport electrons in opposite di- 
rections for opposite spin orientations. At the same time, Luttinger-liquid correlations fractionalize 
electrons injected at a given edge into counter-propagating charge pulses carrying definite fractions 
of the elementary electron charge. 
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A quest for solid-state medium providing means to pro- 
duce, transport, manipulate, and detect quantum entan- 
glement has been fueling tremendous research activity in 
recent years. Many of the practical proposals rely upon 
an innately quantum-mechanical electron spin as an ele- 
mentary building block for quantum bits and ultimately 
quantum computation [1 . The promise for spin-based 
quantum manipulations hinges on its relatively long re- 
laxation and coherence times, in contrast to charge-based 
degrees of freedom that are naturally more susceptible 
to various decoherence processes inherent to a solid-state 
environment. An alternative fault-tolerant charge-based 
quantum computation relics on topological protection of 
braiding statistics of non-Abelian anyons in certain quan- 
tum Hall states [2], although the experimental situation 
there still offers a formidable challenge [3]. In this Letter, 
we propose to utilize spin-dependent topological proper- 
ties of newly-discovered two-dimensional quantum spin 
Hall insulators (QSHI) j3HS] in combination with spin 
entanglement of Cooper pairs (CP's) in superconductors 
as a starting point for spin-based macroscopic quantum 
manipulations [TM12]. Recent experimental strides in re- 
alizing mesoscopic CP splitters [T3J [H] and measuring 
their cross correlations [T5] show the feasibility of our 
proposal. 

While spins appear to be good candidates for initial- 
ization (e.g., by spin-exchange or superconducting corre- 
lations) and local gating (both by magnetic and electric 
fields) of quantum bits, the transport and detection re- 
main more challenging. For the latter, some kind of spin- 
to-charge conversion appears to be necessary at present 
[16]. In this Letter, we utilize topological helicity of the 
QSHI edge states to convert spin entanglement into mea- 
surable charge-current correlations, which are amenable 
to modern experimental capabilities [17] . While we will 
not study here the feasibility of using the same system 
to transport spin entanglement, its potential to this end 
should also be clear from our discussion. Our proposal 
exploits the ideas [T0l - fT2] for injecting CP's into two 
Luttinger-liquid (LL) wires. A key role will be played 



by electron-electron interactions along the edges away 
from the contacts with the superconductor (SC), which 
govern nonchiral LL charge fractionalization and interac- 
tions with Fermi-liquid contacts [18] as well as suppress 
same-edge tunneling [TO] E]. The tunneling cross cor- 
relations thus contain a wealth of information about LL 
effects, SC spin pairing (especially when dealing with un- 
conventional SC's [IS]), and topological properties of the 
QSHI. 

At sufficiently low voltages and temperatures (on the 
scale set by the superconducting gap A), two electrons 
that initially constitute a CP on the superconducting side 
are injected into the QSHI coherently, with a time delay 
of up to ~ A -1 . The larger this delay, the weaker the LL 
effects on the suppression of the tunneling into the same 
edge [TDj. Since this does not affect the low-energy tun- 
neling exponent, however, we can use a simplified model 
of equal-time CP injection [11] . reinserting the delay ef- 
fects [10] via the appropriate energy normalization by A 
when necessary. The (equal-time) tunneling Hamiltonian 
has two pieces corresponding to the CP injection into the 
same (H') or different (fr) edges at x = 0: 



H' = Ye~ i2eVt 
H' = te~ l2eVt 



V^V^W + ^Vf (o) 



H.c., (1) 
H.c., (2) 



where V is the voltage applied between the SC and 
QSHI and u/l label electrons in the upper/lower edges. 
?/'j u ^' [x) are electronic field operators for the helical edge 
modes. Figure [lja) shows a schematic of our setup. 
Since the CP density of states in the SC is a delta func- 
tion at the Fermi level, the effective tunneling Hamil- 
tonian above determines the total current injected from 
the SC. We will focus on the weak-tunneling limit, where 
the overlap between different CP's injected in the QSHI 
is negligible and the spin-singlet entanglement is maxi- 
mized. 

In Eq. 0, we assumed symmetric tunneling into the 
two edges and in the second equation, we accounted 



2 



(a) SC 

qshi m C 

x=0 



(0 



T 



(b) 



e- 




FIG. 1: (Color online) Schematic of the model, (a) Same- 
edge (H') and cross-edge (H') singlet electron-pair injection 
processes from the superconductor into the upper and lower 
QSHI edges, (b) A single right-moving electron injected into 
an edge LL "pumps" a charge of e(l — g)/2 to the left. The net 
fractionalized charges propagating to the right (left) are thus 
respectively given by e± = e(l ± g)/2 [18] . (c) A schematic 
bending of the edges shows a close analogy of our model to 
the proposed entangled CP injection into carbon nanotubes 
[11] . The essential differences stem from the LL correlations 
along the edges vs those at carbon nanotubes' ends and the 
natural spin filtering provided by the edge-state helicity in 
the QSHI. 



for singlet pairing of the CP's. It should be clear that 
any asymmetry in the two tunneling contacts would en- 
hance the relative role of the same-edge tunneling. The 
electron-electron repulsion, on the other hand, can re- 
verse this tendency toward the different-edge tunneling 
[TU1 [TT] . There are at least two factors, however, that 
may hinder this most interesting scenario. First of all, the 
cross-edge tunneling amplitude T is suppressed exponen- 
tially for the edge separations d larger than the supercon- 
ducting coherence length £. Even more problematically, 
CP injection into different edges is suppressed as a power 
law in fcpii (where Uf is the SC's Fermi wave number) 
due to destructive Friedel-like interference for electrons 
tunneling from the SC into different contacts [TO]. This 
issue can be overcome by injecting CP's via an interme- 
diate layer of material with a longer Fermi wavelength 
and proximity-induced SC gap |13) . 

Electron tunneling into different edges is spin sym- 
metric if the full structure is inversion symmetric. The 
latter can be affected by, e.g., applying a local strain 
or gate voltage to one of the edges, thus locally mod- 
ifying its spin-orbit coupling. We will account for this 
tunable asymmetry by a relative angle 6 between spin- 
quantization axes at the two contacts. See Figs.JIJa) and 
|TJb) . Tunneling Hamiltonian ^ then has to be modified 
by a relative spin rotation at, say, the lower edge: 



-> cos(0/2)<4j. ± sm(0/2)vS 



td 



(3) 



In comparison to Ref. [TT], our 6 in Eq. ^ has a sim- 
ilar effect on the edge-current cross correlations as spin 
filtering in carbon nanotubes along axes misaligned by 
8. The key practical difference here is that such spin 



filtering is difficult and has not yet been achieved experi- 
mentally, while our proposal is based on measuring unfil- 
tered charge currents. In other words, an analog of spin 
filtering is already built into helicity of the edge states. 

The effective edge-state Hamiltonian for the QSHI, in- 
cluding the interband and intraband forward-scattering 
processes near the Fermi points is 



k—u.l 



dx 
l/r 



(4) 



are linear 
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combinations of the "spin-s" 
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Here, </> (fe) = ( 
for k — u/l 

bosonic density operators <pT', such that the bosoniza- 
tion identity for the fermionic field operators along the 
upper edge with the spin- up/down electrons moving to 



the right/left is ipi i (x) oc e IS *="' w . Since the relation- 
ship between the spin and orbital chirality reverses at 
the opposite edge, ipi L \x) oc e^ ls ^ l) ^ x \ In our conven- 
tion, the commutation relations are [9^ k \x), (jft*){x')] — 
(i7r/2)sign(x — x'). g in Eq. Q parametrizes the strength 
of the electron-electron interactions (g = 1 for free elec- 
trons and < g < 1 for repulsive interactions). In- 
terband scattering between the right and left movers 
within a given edge leads to nontrivial correlation effects 
(i.e., j / 1), making the QSHI fundamentally distinct 
from an integer quantum Hall insulator with chiral edge 
states [53]. In the special case of equal-strength inter- 
band and intraband scattering, v = vp/g, in terms of 
the bare Fermi velocity vp. In the general case with 
broken Galilean invariance, however, g and v are inde- 
pendent phcnomenological parameters governed by the 
interplay between band-structure and correlation effects 
[2D] , A rough estimate [55] gives g > 0.5 for helical edge 
states in a HgTe quantum well [5] , corresponding to mod- 
erate interaction effects while the two-electron disorder 
and umklapp backscattering are still irrelevant |21j . 

In order to measure current cross correlations, we place 
current sensors along the edges at points 1, 2, 3, and 4, as 
shown in Fig.JTJa). Alternatively, we can divert edge cur- 
rents into voltage probes [3J and measure low-frequency 
voltage correlations (SViSVj). Note, however, that the 
Fermi-liquid voltage probes would affect incoming cur- 
rents, effectively masking the LL charge fractionalization 
at low voltages 18J. Using bosonic representation for the 
current operators along the edges, I 1 -^ — vgd x 9^ k ' /it, we 
evaluate the symmetrized noise spectrum 



Sij(uj) = Sji(-uj) = 



dte™ ({Sliit), SljiO)}) , (5) 



where 8Ii(t) = — (Ii(t)). i labels four outgoing chan- 
nels: i = 1, upper right; 2, upper left; 3, lower left; and 
4, lower right branches. To this end, we use standard 
Keldysh formalism [5T] . Note that when calculating noise 
correlations to the leading (second) order in the tunneling 



3 



coefficients 1^, we do not distinguish between the noise 
(SliSIj) and current (Iilj) correlators, since the differ- 
ence enters at the fourth order in tunneling. The 
current direction is always chosen away from the SC. 

Let us first discuss the most interesting regime, when 
a CP injected from a SC tunnels into two separate edges. 
This is governed by Eq. ([2| and sketched by the cross- 
edge right- and left-moving electron pairs in Fig.jlja). A 
finite 9, furthermore, leads to correlations between leads 
1 and 3 as well as 2 and 4. Initial spin-singlet entan- 
glement is thus converted into the ^-dependent cross- 
edge current correlations. In addition to this, there is 
a purely LL entanglement stemming from the injected 
charge fractionahzation [18] . which has to be included 
on equal footing with the spin entanglement. This is 
sketched in Fig. [ljc): A single electron injected into a 
given edge sets off a counterpropagating nondispersive 
charge wave, which eats up a fraction (1 — g)/2 of the 
elementary electron charge e. 

Putting these superconducting and LL correlations to- 
gether, we find for the low-frequency noise cross correla- 
tors Si4 = 5*23 = S+: 
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We are focusing on the low-temperature regime, fc^T <C 
eV, so that thermal noise can be neglected, ep is the 
characteristic bandwidth of the LL, typically of the order 
of the Fermi energy, and 7 = (g+g^ 1 ) /2— 1 is the spinless 
nonchiral LL bulk exponent for a single-electron injec- 
tion [5D]. Apart from the generic LL tunneling anomaly, 
there are two types of interesting factors in Eq. Q: LL 



factors el 



e 2 (l + <7 2 )/2 and 2e + e_ =e 2 (l- 5 2 )/2 



reflecting charge fractionahzation [18] and trigonometric 
factors cos 2 (#/2) and sin 2 (6>/2) reflecting spin-singlet en- 
tanglement. Similarly, we find for S13 — S24 = S_: 
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(e 2 + e 2 ) sin 2 - + 2e+e_ cos 2 
2 1 — g 2 cos 9 



(7) 



For the total interedge correlation S u i = (I u Ii) — Si 3 + 
S14 + S 23 + S 24 = 2(S+ + 5_), we have S ut /V oc V 2 ~<, 
independent of the charge fractionahzation and angle 9. 
The self-correlators are determined by the respective av- 
erage currents, as generally expected for the Poissonian 
statistics in the tunneling regime: 



(e\ + e 2 _){h)/e 



fl 



(8) 



In the special case when 9 = and g = 1, there is a 
perfect positive cross correlation of currents at contacts 



1, 4 and 2, 3 (i.e., S+ = So) and zero correlation at 
contacts 1, 3 and 2, 4 (i.e., S_ =0). The former is cer- 
tainly a manifestation of many-body entanglement, since 
nonintcracting electrons are required to have nonposi- 
tive cross correlations in general multiterminal structures 
[17] , According to Eqs. ^ and |7|, the LL parameter 
g can be measured in the inversion-symmetric configura- 
tion (so that 9 = 0) using 



9 = 



= ^/(s + -s_)/(s+ + s_). 



(9) 



If 9 is unknown and not precisely controlled, but can, 
nonetheless, be swept over the half-period (0,tt), g will 
be found by maximizing the absolute value of the right- 
hand side in Eq. which equals g\J cos 9 according to 
Eqs. fl6J, ([T]). The same g experimentally extracted from 
Eqs. Q and ^ would provide direct evidence of spin 
entanglement as manifested by the interedge correlations 
(even if g turns out to be trivial, i.e., close to 1). 

Detrimental to the above formulation would be any 
backscattering along the edges. While the time-reversal 
symmetry protects against backscattering on nonmag- 
netic impurities, there would inevitably be interactions 
of the edge modes with nearby Fermi liquids (such 
as electrostatic gates, metallic contacts etc.), and the 
edges would terminate somewhere, possibly merging with 
a Fermi-liquid continuum. As was already mentioned 
above, the Fermi-liquid probes measuring current fluc- 
tuations would conceal the LL fractionahzation [TS], ef- 
fectively removing positive cross correlations associated 
with the LL physics. In this case, we would expect g 
in Eqs. (|6])-([9| to be close to unity. More problemati- 
cally, backscattering engendered by charge fluctuations 
on nearby gates or other metallic regions, which would 
scramble ballistic propagation of edge states, needs to be 
eliminated from a practical implementation of our pro- 
posal. It should be clear, however, that various backscat- 
tering and boundary-related effects should not affect the 
finite-frequency correlations on time scales shorter than 
those at which they set in while still longer than A -1 , so 
that the CP tunneling delay is not resolved. 

In the absence of the LL fractionahzation, g = 1, the 
noise S± oc (1 ± cos#)/2 has the same angular depen- 
dence as the spin-current noise in the carbon-nanotube 
proposal of Ref. [11 . The 9 dependence here is the same, 
in turn, as in the Einstein-Podolsky-Rosen thought ex- 
periment [25] on a delocalized spin-singlet state of two 
electrons. It is instructive to understand a close relation- 
ship between our setup and that of Ref. [TT] by schemat- 
ically folding two sides of cither infinite QSHI edge into a 
semi-infinite "wire," as sketched in Fig. [TJb). This wire is 
then formally identical to a semi- infinite spinful LL, only 
without interactions between spins up (down) moving in 
opposite directions as well as between spins up and down 
moving in the same direction. There is still, however, a 
nontrivial interaction between spins up and down mov- 



4 



ing in opposite directions, which corresponds to g =/= 1 in 
our original model shown in Fig. [ija) . Hypothetical spin 
filtration pj] in semi-infinite carbon nanotubes is thus 
provided in our system by effectively unfolding the ID 
modes of such semi-infinite wires into infinite edges, with 
half the number of ID bands and helical spin character 
(i.e., spins up and down going in opposite directions). 

Owing to the difference in the interaction channels here 
and in Ref . , we find a hierarchy of the tunneling ex- 
ponents that is more similar to that discussed in Ref. [TO] 
for infinite LL wires. Injecting a CP into opposite QSHI 
edges is characterized by the exponent 27 = g + g^ 1 — 2, 
as has already been discussed above [Eqs. ^ and Q]. 
We find that the same-edge tunneling of a singlet CP is 
suppressed by an additional factor ~ (eVyA) 2 ^ 7 ~ 7 ) with 
respect to Eqs. |6| and ([7]), where 7' = g^ 1 — 1 is the end- 
tunneling exponent for a nonchiral spinless LL |20j . The 
fact that 7' > 7 for g < 1 insures that the split-tunneling 
dominates at low energies P33IH] over the same-edge tun- 
neling, which is parasitic for our purposes. In particular, 
the same-edge tunneling would introduce large positive 
cross correlations S12 and S34, while not contributing to 
the interedge correlations. 

The situation is quite different in the case of unconven- 
tional SC's with triplet spin pairing 19]. For example, 
a CP in the equal-spin pairing phase could be injected 
into the same (left- or right-moving) edge mode, in the 
case of the same-edge tunneling, which is governed by the 
LL exponent 47 + 2 vs 27 for the cross-edge tunneling. 
The relative suppression factor here is ~ (e^/A) 27+2 , 
with the quadratic factor stemming from the Pauli ex- 
clusion that hinders injection of two same-spin electrons 
into the same edge band [26]. For a finite 8, there are also 
nontrivial interedge correlations, dictated by the super- 
conducting spin pairing. In this context, we can specu- 
late, furthermore, that a QSHI not only provides a fruit- 
ful medium for injecting and subsequently manipulating 
spin-entangled CP's but can also serve as a probe for spin 
pairing, which can thus be manifested via current cross 
correlations along the edges. This offers an alternative 
to the proposal of experimentally probing the supercon- 
ducting order parameter via nonlocal spin and charge 
pumping [27j . 

We finally note that our setup is distinct from that of 
the Hanbury Brown and Twiss correlations of CP's in 
the QSHI discussed recently in Ref. 28 . In the latter 
implementation, CP's are injected ballistically from an 
edge portion that opens a gap due to a strong proxim- 
ity effect to a nearby SC. The induced superconductivity 
provides perfect Andreev reflection for electrons incom- 
ing from the normal region of the edge. Time-reversal 
symmetry prohibits normal scattering in this situation. 
This Andreev reflection is found to induce ordinary (i.e., 
negative) current cross correlations with the second edge 
that is connected to the first by a point contact [25] ■ In 
the present situation, on the other hand, the contact be- 



tween the SC and the QSHI is weak and the electrons 
are thus passing along the edges almost ballistically with 
little Andreev scattering. Those electrons that rarely un- 
dergo Andreev reflection, in turn, scatter to a different 
edge as holes, which is called crossed Andreev reflection 
(while the same-edge Andreev reflection is suppressed by 
the LL correlations). This effective flipping of the scat- 
tered carrier charge reconciles our positive cross correla- 
tions with the Hanbury Brown and Twiss perspective of 
Ref. [25]. 
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